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ABSTRACT 



Context. Evidence of magnetic interaction between late-type stars and close-in giant planets is provided by the observations of stellar hot spots 
rotating synchronously with the planets and showing an enhancement of chromospheric and X-ray fluxes. Possible photospheric signatures of 
such an interaction have also been reported. 

Aims. We investigate star-planet interaction in the framework of a magnetic field model of a stellar corona, considering the interaction between 
the coronal field and that of a planetary magnetosphere moving through the corona. This is motivated, among others, by the difficulty of 
accounting for the energy budgets of the interaction phenomena with previous models. 

Methods. A linear force-free model is applied to describe the coronal field and study the evolution of its total magnetic energy and relative 
helicity according to the boundary conditions at the stellar surface and the effects related to the planetary motion through the corona. 
Results. The energy budget of the star-planet interaction is discussed assuming that the planet may trigger a release of the energy of the coronal 
field by decreasing its relative helicity. The observed intermittent character of the star-planet interaction is explained by a topological change of 
the stellar coronal field, induced by a variation of its relative helicity. The model predicts the formation of many prominence-like structures in 
the case of highly active stars owing to the accumulation of matter evaporated from the planet inside an azimuthal flux rope in the outer corona. 
Moreover, the model can explain why stars accompanied by close-in planets have a higher X-ray luminosity than those with distant planets. 
It predicts that the best conditions to detect radio emission from the exoplanets and their host stars are achieved when the field topology is 
characterized by field lines connected to the surface of the star, leading to a chromospheric hot spot rotating synchronously with the planet. 
Conclusions. The main predictions of the model can be verified with present observational techniques, by a simultaneous monitoring of the 
chromospheric flux and X-ray (or radio) emission, and spectropolarimetric observations of the photospheric magnetic fields. 
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1. Introduction 

More than 350 extrasolar giant planets are presently knowrfl, 
among which ~ 25 percent have a projected orbital semi- 
major axis smaller than 0.1 AU. Such planets are expected 
to interact with their host stars, not only through tides, but 
also with other mechanisms, possibly associated with magnetic 
fields. Specifically, reconnection between the coronal field of 
the host star and the magnetic field of the planet is expected 
to release heat, produce hydromagnetic waves, and acceler- 
ated particles that may be conveyed onto the stel lar chromo- 
sphere produc i ng a localized enhanced emission dC untz et alj 
2000 !: 'ip et al.l l2004, I Preusse etal]l20 06l: 'Mc lvor et aL,200a 
Cranmer & Saar 2007). Indeed, Ishkolnik et al.l tOOj BOOsf) 
show that HD 179949 and v Andromedae have chromospheric 
hot spots that rotate with the orbital periods of their inner plan- 
ets. Two other stars, HD 189733 and r Bootis, show evidence of 
an excess of chromospheric variability, probably due to flaring, 
that is modulated with the orbital periods of their respective 
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planets (cf. also IWalker et al Modelling such features 

can allow us to obtain information on the magnetic fields of th e 
planets, although in an indirect way (cf. IShkolnik et a"l]|2008h . 
which is of great importance to constrain their internal struc- 
ture and evolution as well as to characterize a possible habit- 
ability in the case of rocky planets. A recent step forward in 
the understanding of the star- planet magneti c interaction (here- 
inafter SPMI) is the work of ILanzal (l2008b . It applies force- 
free and non-force-free coronal field models to account for the 
observed phase lags between planets and synchronous chro- 
mospheric hot spots. Moreover, it presents conjectures about 
a possible influence of a close-in planet on the hydromagnetic 
dynamo action occurring in its host star. In the present paper, 
we shall extend that work discussing topology, total energy, 
and relative magnetic helicity of the coronal field and inves- 
tigating their role in the effects associated with the motion of 
a hot Jupiter inside a stellar corona. We shall propose a new 
mechanism to account for the energy budget of SPMI and its 
intermittent nature, as revealed by the latest observations (see 
Sect.|2]i. Moreover, we shall discuss how the different coronal 
field topologies and the presence of a planet embedded in the 
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corona may affect the formation of prominence-like structures 
as well as the X-ray and radio emissions of a star. 

2. Observations 



Evidence o f SPMI was reported by IShkolnik et alJ (12003 
20051 l2008h who observed chromospheric hot spots rotating 
with the orbital period of the planets in HD 179949 and v And 
instead of their respective rotation periods. The spots were not 
located at the subplanetary point, i.e., along the line joining the 
centre of a star with its planet, but lead the planet by ~ 70° 
in the case of HD 179949 and ~ 170° in the case of v And, 
respectively. The flux irradiated by the hot spot was greater in 
the case of HD 179949, with an excess power of ~ 10^" W, cor- 
responding to that of a large solar flare. The hot spot synchro- 
nized with the planet was observed in four out of six seasons 
in HD 179949, suggesting that it is not a steady phenomenon, 
but may have a lifetime not exceeding 300 - 400 days. A sim- 
ilarly intermittent interaction is suggeste d by the observations 
of V And, leading IShkolnik et alj (120081) to propose that these 



on/off SPMI transitions are related to a variation of the con- 
figuration of the coronal magnetic field of the host stars, pos- 
sibly connected to their acti vity cycles. For the p lanet hosts 



HD 189733 and HD 73256. IShkolnik et al.1 (12008*) suggest a 



correlation of the amplitude of the intranight variability of the 
Ca II K line core flux with the orbital phase of the planet, with 
a maximum of activity leading the planet by ^ 70° in the case 
of HD 189733. This may be due to some chromospheric flaring 
activity synchronized with the planet. 

The X-ray flux coming from HD 179949 appears to be 
modulated with the o rbital period of t he planet with a varia- 



tion of » 30 percent dSaar et alj 120070 . A statistical analysis 



of a sample of late-type stars reveals that those hosting a hot 
Jupiter closer than 0.15 AU have, on the average, an X-ray flux 
ss 3 - 4 times greater than stars hosting planets farther than 1 .5 
AU dKashvap et alj|2008l) . This suggests that a close-in giant 
planet may increase the X-ray activity level of its host star 

Circumstantial evidence of planet-induced photospheric 
magnetic activity in s tars hosting hot Jupiters has been re- 
ported in T Bootis bv IWalker et all (12008.) and in CoRoT-2 a 
and CoRoT-4a by iLanza et alj (l2009ai nPagano et all (l2009l) . 
and Lanza et al. I (l2009bl ). respectively. For r Boo there is evi- 
dence of an active region leading the planet by about 70° whose 
signatures have been observed both in the Ca II K line core flux 
and in the wide-band optical flux as monitored by the MOST 
(Microvariability and Oscillation of STars) satellite. The region 
shows a rotational modulation of ~ 1 mmag. In 2004 it re- 
sembled a dark spot of variable depth, while in 2005 it varied 
between bright and dark. Since t Boo has an average rotation 
period of 3.3 days, synchronized with the orbital period of the 
planet, and displa ys a surfac e differential rotation comparable 

to that of the Sun dCatala et al 2007) . it is the constancy of the 

phase shift between the spot and the planet that gives support 
to a possible SPMI in this case. As a matter of fact, the hot spot 
can be traced back to 2001, thanks to a previous Ca II K line 
flux monitoring. 

The planet host CoRoT-4a as a spectral type similar to t 
Boo and its rotation appears to be synchronized on the av- 



erage with that of its planet. Modelling two month s of unin- 
terrupted CoRoT observations, iLanza et al. I (l2009bl) found ev- 
idence of a persistent active region at the subplanetary lon- 
gitude. Five months of observations of CoRoT-2a revealed a 
modulation of the total spotted area with a cycle close to ten 
synodic period s of the hot Jupiter with respect to the stellar ro- 
tation period (Lanza et alj|200 9a)). again suggesting some kind 
of SPMI. Pagano et al. ( 2009) found that the variance of the 
stellar flux was modulated in phase with the planetary orbit, 
with minimum variance at phase 0.3 - 0.4 and maximum at 
phase 0.8 - 0.9. Finallv. iHenrv et alj (120021) found evidence of 



a photospheric spot in HD 192263 that in two observing sea- 
sons rotated with th e orbital period of it s planet for at least three 
rotation cycles (see Santos et"al]|2003 ). 



This preliminary evidence of photospheric cool spots syn- 
chronized with a hot Jupiter is impossible to explain in the 
framework of a magnetic reconnection model because recon- 
nection leads to energy release and thus heating of the atmo- 
sphere. Conversely, it suggests that the planet may affect in 
some way the stellar dynam o action or th e emergence of mag- 
netic flux, as conjectured by lLanzal ( 2008 ). 

Spectropolarimetric observations can be applied to map 
photospheric fields outside starspots and indeed for t Boo a 
sequence of maps has been obtained that suggests that the mag- 
netic activity cycle of th e star is as short as ~ 2 years, instea d 
of 22 years as in the Sun dDonati et al.l2008l:lFares et al.ll2009l) . 

Solar system giant planets emit in the radio domain, mainly 
via the electron-cyclotron maser mechanism. Most of the 
power is coherently radiated near the cyclotron frequency /^e = 
2.8B MHz, where B is the intensity of the planet's magnetic 
field in Gauss. In the case of Jupiter, the maximum of the flux 
falls around 40 MHz. Solar system planets follow a scaling 
law that relates their emitted radio power to the power supplied 
by the impinging solar wind at their magnetospheric boundary. 
Generalizing that scaling law, it is possi ble to predict exoplan- 
etary radio emission powers (cf., e.g. , IStevensI 120051: IZarka 
2007 : Jardine & Collier CameronI 2008 ). For r Boo and some 



other systems, fluxes between 30 and 300 mJy, i.e., within the 
detection limits of some of the largest radio telescopes, have 
been predicted for favourable conditions. Nevertheless, no pos- 
itive detection has been reported yet, not only in the case of t 
Boo, observed at several epochs at 74 MHz wit h upper limits 
between 135 and 300 mJy ( Lazio & FarrelJ2007 ), but also for e 
Eridani and HD 12831 1 at 150 MHz (' George & Stevensll2007l) 
for which tight upper limits of 10-20 mJy were derived, and for 
HD 189733, for which an upper Umit of ~ 80 mJy was reached 



over the 307-347 MHz range (ISmith et al.1 120091) . It is hoped 
that the next generation of low-frequency radio telescopes can 
lower those limits by at least one order of magnitude allowing 
us to clarify whether the missed detections are due to a lack of 
sensitivity or a lack of emission from the exoplanets at those 
frequencies. As a matter of fact, if hot Jupiters have magnetic 
field strengths comparable to that of Jupiter, i.e., B ~ 14.5 G 
at the poles (with B ~ 4.3 G at the planet's equator), most of 
their radio emission falls at frequencies lower than those of the 
above mentioned surveys. Another possibility is that the emis- 
sion is beamed out of the Une of sight, or that the stellar wind 
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power (kinetic and/or magnetic) is significantly smaller than 
assumed. 



3. Coronal field model 

We adopt a spherical polar coordinate frame having its origin at 
the baricentre of the host star and the polar axis along the stellar 
rotation axis. The radial distance from the origin is indicated 
with r, the colatitude measured from the North pole with 0, and 
the azimuthal angle with (p. The planet orbit is assumed circular 
with a semimajor axis a and lying in the equatorial plane of the 
staid 

Close to a star, the magnetic pressure in the corona is much 
greater than the plasma pressure and the gravitational force, 
so we can assume that the corona is in a force-free magneto- 
hydrostatic balance, i.e., the current density / is everywhere 
parallel to the magnetic field B, viz. J x B - 0. This means 
that V X B = aB, with the for ce-free parameter a constant 
along each field line ('P riestlll982h . If a is uniform in the stellar 
corona, the field is called a linear force-free field and it satisfies 
the vector Helmoltz equation V^B + a^B = 0. Its solutions in 
spheri cal geometry have been studied by, e. g., lchandrasekhai 
( Il956l) andlChan drasekhar & Kendall ( Il957l) . 

Linear force-free fields are particularly attractive in view 
of their mathematical symplicity and their minimum-energy 
properties in a finite domain, as shown by Woltjer (1958). 
Specifically, in ideal magnetohydrodynamics, the minimum en- 
ergy state of a magnetic field in a finite domain is a linear 
force-free state set according to the boundary conditions and 
the constrain posed by the conservation of the magnetic he- 
licity. In a confined stellar corona, magnetic dissipation is lo- 
calized within thin current sheets whose global effect is that 
of driving the field configuration toward the minimum energy 
state compatible with the conservation of the total helicity, the 
so-c alled Taylor's state, which is a linear fo rce-free field (see, 
e.g., iHevvaerts & Priestlll 984 iBergedl 19851 for details). Only 
in large flares the total helicity is not conserved due to strong 
turbulent dissipation and ejection of magnetized plasma. 

The total helicity of a confined magnetic structure con- 
strains its free energy, that is the energy that can be released in 
a magnetic dissipation process. It is the difference between the 
initial energy of the field, which is usually in a non-linear force- 
free state (i.e., with a non-uniform a), and the energy of the 
linear force-free field satisfying the same bo undary conditions 
and h aving the same total helicity (see, e.g., iRegnier & Priest 
20071). Therefore, a decrease of the total helicity produced by a 
change of the boundary conditions or a strong magnetic dissi- 
pation will in general make available more free energy for the 
heating of the corona. 

To model the magnetic interaction between the stellar coro- 
nal field and a close-in planet, we consider only the dipole-like 
component (i.e., with a radial order n - 1) o f the l inear force- 
free solution of £handrasekhar & Kendalll ( Il957b because it 
has the slowest decay with distance from the star and therefore 



^ About 70 percent of the hot Jupiters within 0. 1 AU from their host 
stars have a measured orbital eccentricity lower than 0.05, consistent 
with a circular orbit (see http://exoplanet.eu/ 1. 



leads to the strongest interaction. Moreover, since the observa- 
tions of star-planet interaction show a hot spot rotating with the 
orbital period of the planet and do not show the periodicity of 
stellar rotation, an axisymmetric field (i.e., with an azimuthal 
degree m = 0) is a good approximation to describe the interac- 
tion, as discussed by lLanzal (l2008l) . 

Our li near force- f ree fie ld can be expressed in the for- 
mulism of Fiver et al. (2004) as: 



B = 



1 



rsinO 



-— r - —0 + aAd> 

r 89 dr 



(1) 



where A{r, 6) is the flux function of the field. Magnetic field 
lines lie over surfaces of constant A(r, 9), as can be deduced 
by noting that B ■ VA = 0. The flux function for our dipole- 
like field geometry is A(r, 9) - BoR^giq) sin^ 9, where 2Bo is 
the magnetic field intensity at the North pole of the star, R the 
star's radius and the function g{q) is defined by: 

[boJ-3/2(q) + coJ3/2(q)] y/q 



[boJ-iiiiqi)) + co73/2(g'o)] V*) ' 



(2) 



where bo and cq are free constants, 7-3/2 and 73/2 are Bessel 
functions of the first kind of order -3/2 and 3/2, respectively, 
q = \a\r, and ^0 = \a\R. Making use of Eq. ([TJ, the magnetic 
field components are: 

R^ 

Br = 2Bo—g{q)cos9, 

Be = -BQ\a\—g'{q)sm9, (3) 
r 

R^ 

Bs = aBo — g{q)sm9, 
r 

where g'{q) = dg/dq. A linear force-free field as given by 
Eqs. (O extends to infinity with an infinity energy. We consider 
its restriction to the radial domain qo < q < q^, where qi^ is the 
first zero of g{q), in order to model the inner part of the stel- 
lar corona where mag netic field lines are cl osed, as discussed 
in Lanzal ( 2008,) (see Chandrasekhai 19561 for the boundary 
conditions at r = tl = qh/la])- 

The magnetic field geometry specified by Eqs. ([3]) depends 
on two independent parameters, i.e., a and bo/cQ. They can 
be derived from the boundary conditions at the stellar photo- 
sphere, i.e., knowing the magnetic field B^''\8,(/)) on the sur- 
face at r = R. Using the ortho gonahty properties of the basic 
poloidal and toroidal fields (see Chandrasekhati 1961 ), we find: 



871 , 

-jBoR^ 



a\BoR'g'iqo) 
—aBoR^ 



j 
j 



cos 9dl., 



B^g^ sin 9d'L, 



B^'^ sin 0^/2, 



(4) 



where S(/?) is the spherical surface of radius R, and dl, - 
R^ sin 9d9d4>. 

The magnetic energy E of the field confined between the 
spherical surfaces r - R and r - ri^ can be found from Eq. (79) 
in § 40 of Chandrasekhar (1961): 



£ = £p {2 + qoqL[8'(qL)f - ql[g'(qo)f - ql} , 



(5) 
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where Ep = is the energy of the potential dipole field 

with the same radial component at the surface r - R, and fi is 
the magnetic permeab ility. The relative magnetic helicity //r, 
as defined by Bergerl (fl985 ). can be found from his Eq. (19) 
and is: 



By 



2g\qo) + 



Stt E 



3 qoEp 



Note that the field obtained by changing the sign of a has the 
same poloidal components and Bg, and energy E, while the 
toroidal component B^ and the relative helicity become op- 
posite. 

In the limit q ^ I, i.e., close to the stellar surface, neglect- 
ing terms of the order O(q^), we find that g{q) ~ qo/q and 

g'iq) qo/q^, independently of bo and cq. From Eqs. (|3]l, we 

see that and Bg close to the star are similar to the analogous 
components of the potential dipole field with the same radial 
component at the surface. 

It is also interesting to study the limit a — > of our force- 
free field. In this limit tl — > oo and the radial extension of the 
field grows without bound. Nevertheless, the field tends to be- 
come indistinguishable from the coiTesponding potential field 
and E — > Ep. In the same limit, the relative magnetic helicity 
grows without bound because — > in Eq. but the he- 
licity densit y (i.e., the magnetic helicity per uni t volume) tends 
to zero (see lBergeilll985[ IZhang & Lowll2005l for further de- 
scription of this limit state). 

For a finite a, E > Ep because the potential field has the 

(s) 

minimum energy for a given BY'. If we consider all magnetic 
fields with one end of their field lines anchored at r = /? and 
the other out to infinity, satisfying the same boundary condi- 
tions of our field at r = R, the field with the lowest possible 
energy is called t he Aly field and its energy ^Aiy - 1 •66£'p (see 



Flyer et al.ll2004h . We assume that the Aly energy is an upper 
bound for the energy of our field because it is the lowest energy 
allowing the field to open up all its lines of force out to infinity 
driving a plasma outflow similar to a solar coronal mass ejec- 
tion. 

From a topological point of view, the fields obtained from 
Eqs. (|3]l can be classified into two classes. If the function g{q) 
decreases monotonously in the interval qo < q < qh, all field 
lines are anchored at both ends on the boundary r - R. On 
the other hand, if the function has a relative minimum (and 
a relative maximum) in that interval, the field contains an az- 
imuthal rope of flux located entirely in the r > R space and 
running around the axis o f symmetry. T he magnetic field con- 
figurations considered by iLanzal ( 120081) to model SPMI are of 
the first kind; an example of a field configuration containing 
an azimuthal flux rope will be discussed in Sect. 14.21 An anal- 
ogous topological classification hol ds in the case of th e non- 
linear force-free fields considered by lFlyer et al. (2004). 

Note that the photospheric magnetic field components can 
be measured by means of spectropolarimetric techniques if the 
star rotates fast enough (v sin / > 10 - 15 km s~') as shown in 
the case of, e.g., r Boo by Catala et alJ (120071) and lDonati et al 



( 2008h . Therefore, Eqs. (@| can be applied to derive the param- 
eters of the coronal field model and its topology, as we shall 
show in Sect. 14.21 



4. Applications 

In this Section we shall consider some applications of the above 
model for the stellar coronal field to the phenomena introduced 
in Sect.|2] 



(6) 4. 1 . Power dissipated in SPMI 



The power needed to explain the excess flux from a chromo- 
spheric hot spot or an X-ray emission synchronous with the 
planet is of the order of 10^" - 10^' W. Magnetic reconnection 
between the stellar coronal field and the planetary field at the 
boundary of the planetary magnetosphere is not sufficient to 
account for such a power, as we show in Sect. l4.1.TI Therefore, 
we propose an interaction mechanism that may be capable of 
sustaining that level of power in Sect. 14.1.21 

4.1 .1 . A simple magnetic reconnection model 

The first mechanism proposed to account for the energy budget 
of SPMI is reconnection between the planetary and the stel- 
lar magnetic fields at the boundary of the planetary magneto- 
sphere. This boundary is characterized by a balance between 
the magnetic pressure of the coronal field and that of the plan- 
etary field. The ram pressure is negligible because the planet 
is inside t he region where the stellar wind speed is subalfvenic 



(cf., e.g., iPreusse et al.l 120051) . and the orbital velocity of the 



planet is about one ord er of magnitude smaller than the Alfven 
velocity (lLanzall2008b . Therefore, assuming a planetary field 
with a dipole geometry, the radius of the planetary magneto- 
sphere R„i, measured from the centre of the planet, is given by: 



Rm — Rn 



B(a,^) 



B, 



pi 



(7) 



where Rp\ is the radius of the planet, B(a, |) is the coronal field 
of the star on the stellar equatorial plane at r = a, and Bpi is the 
magnetic field strength at the poles of the planet. 

We specialize our considerations for the SPMI model of 
HD 179949 developed by Lanza (2008) that assumes bo/c^ = 
-1.1 and a - -0.12/? ' to explain the phase lag between the 
chromospheric hot spot and the planet. The magnetic field in- 
tensity vs. the distance from the star is plotted in Fig. [T| The 
field strength decreases as that of a potential dipole field close 
to the star, i.e., where g{q) ~ qo/q (cf. Sect.O, then it decreases 
more slowly for r > (6 - 1)R because the decrease of g{q) be- 
comes less steep (cf., e.g., Fig.|3]i. 

At the distance of the planet, i.e., a - 1J2R, the field is 
reduced by a factor of ~ 300 with respect to its value at the 
stell ar surface. Assumin g a mean field at the surface Bo = 10 G 
(cf. iDonati et al ] l2008h . the field at the boundary of the plane- 
tary magnetosphere is B( a, ^) ~ 0.03 G. Adopting a planetary 
field Bpi = 5 G (cf., e.g., iGrieBmeier eT al.ll2004 and a radius 
of the planet equal to that of Jupiter, we find Rj„ ~ 5.5/?pi - 
3.8 X 10** m. The power dissipated by magnetic reconnection 
can be estimated as: 



Pa^7-[B(a,-)rRiv,,u 

fi L 



(8) 
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where < y < 1 is a factor that depends on the angle b etween 



Priest! 12003). 



the interacting magnetic field lines (see, e.g., 
and Vrei is the relative velocity between the planet and the stel- 
lar coronal field. Adopting the parameters for HD 179949 and 
y = 0.5, we find Fd - 13 x 10'^ W, which is insufficient by 
at least a factor of 10^ Since Pd ~ B'*^■'Bp,''^ to explain the 
observed power the stellar surface field should be increased to 
Bq ~ 180 G, which is too hig h in view of the me an values mea- 



sured at the surface of r Boo ( Donati et alj |2008) that is a faster 
rotator than HD 179949, but with a comparable X-r ay lumi- 
nosity dKashvap et al.ll2008h . lChristensen et al.l (l2009h suggest 
that the surface fields of hot Jupiters may be up to 5 - 10 times 
stronger than that of Jupiter. However, if we adopt Bpi - 75 G 
and consider the dependence, we obtain a power increase 
by only a factor of 6. 

4.1 .2. Magnetic energy release induced by a liot 
Jupiter 

To solve the energy problem, it is important to note that the en- 
ergy released at the reconnection site is only a fraction of the 
energy that is actually available. The emergence of magnetic 
flux from the convection zone and the photospheric motions 
lead to a continuous accumulation of energy in the coronal 
field, independently of the presence of a planet. Our conjec- 
ture is that a close-in planet triggers a release of such an ac- 
cumulated energy, modulating the chromospheric and coronal 
heating with the orbital period of the planet. 

As a matter of fact, the reconnection between the coronal 
and the planetary field lines has remarkable consequences for 
the topological structure of the coronal field. The final con- 
figuration of the coronal field after reconnection is in general 
a non-linear force-free one, i.e., a is no longer spatially uni- 
form. This happens because the field must simultaneously sat- 
isfy the boundary conditions at the photosphere, which impose 
the value of a close to the star, and those at the reconnection 
site close to the planet which in general will not be compati- 
ble with a uniform value of a. Since a must be constant along 
a given field line, the field re-arranges itself into a configura- 
tion with different values of a along different field lines to sat- 
isfy the boundary conditions. We know from Woltjer theorem 
that the energy of the field in such a state is greater than that 
of the linear force-free field with the same magnetic helicity 
and boundary conditions at the stellar surface. In other words, 
some energy can be released if the field makes a transition to 
this linear force-free state, restoring its unperturbed configura- 
tion. Such a transition can be stimulated by the fact that the 
helicity in the region where the field has reconnected is lower 
than the initial helicity because reconnection processes lead to 
a steady dissipation of helicity. Considering the volume V{t) 
occupied at any given time t by the plasma where reconnection 
has just occurred, the rate of helici t y chan ge inside V(t) is given 
by Eq. (13) of He vvaerts & Priesl (Il984 : 



Dt 



{A-v,,i)(B-dS)+ I ^^^dS- 

S(V) Js(V) 



J ^-dS-lf- 

Jsm o" Jy 



-dV, (9) 



where S(V) is the surface bounding the volume V(t), A is the 
vector potential of the magnetic field, viz. V x A = B, and cr 



is the electric conductivity of the plasma. We assume that the 
reconnected field is in a force-free state, thus - aB, and 
choose a gauge transformation for the vector potential to have 
A = a^^B. Since Vrei ■ B ~ in the reconnected region, we 
have: 



r 

—dV, 
Jv cr 



(10) 



where (ff"' ) is some average value of a^^ over the volume V{t). 
Eq. ( fTOl ) shows that the absolute value of the helicity decreases 
steadily inside the volume V{ f) following the motion of the re- 
connected region through the stellar corona. 

We conjecture that once a magnetic energy release is trig- 
gered by such a decrease of helicity in the reconnection vol- 
ume, it extends to the whole flux tube connecting the planetary 
magnetosphere with the stellar surface and proceeds faster and 
faster thanks to a positive feed-back between energy release 
and helicity dissipation. In other words, the fast and localized 
energy release produces a turbulent plasma which in turn en- 
hances turbulent dissipation of magnetic energy and helicity 
giving rise to a self-sustained process. 

An order-of-magnitude estimate of the helicity dissipation 
rate can be based on Eq. ( fTOb . now considering turbulent dissi- 
pation; it can be recast in the approximate form: 



DIh. 
Dt 



(11) 



where 6 is the lengthscale of variation of the magnetic field 
associated with the current density, i.e., J ~ Blidfi), and 
is the dissipated power. We assume that the largest turbulent 
lengthscale is (5 ~ \a\-^ = 0.127?-' = 6.7x 10'' m forT? = 8x10** 
m. Adopting the maximum dissipated power suggested by the 
observations, i.e., f d = 10^' W, we find: 



DIh. 
Dt 



8.37 X lO^'* T^ m" s" 



(12) 



In the above derivation we have guessed the dissipated 
power from the available observations. For the mechanism to 
be at least plausible, we need to show that it is possible to 
sustain that level of dissipated power. This implies a compu- 
tation of the energy stored in the non-linear force- free field that 
is produced by the reconnection process, which is beyond the 
capability of our simple linear model. Nevertheless, we may 
estimate the energy available by means of a simplified argu- 
ment that is by no means rigorous, but has the advantage of 
using our model for linear force-free configurations. Its results 
should be regarded only as an illustration of the plausibility of 
the proposed mechanism, deferring a more detailed and rigor- 
ous treatment to future studies. 

To compute the available energy, we assume that the recon- 
nection between the coronal and planetary fields produces a 
dissipation of the helicity of the initial field configuration, i.e., 
that unperturbed by the planet. Again, we specify our consid- 
erations for th e case of HP 179949, adopting the linear force- 
free model of iLanzal (l2008h . Therefore, we fix = -0.127? ' 
and boico - -1.1 in the initially unperturbed state, giving a 
total energy EjE^ - 1.0273 and a relative helicity |//r| - 
55.3116 Byi'^ (cf. Eqs.|5]and|6]l. The boundary conditions at 
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the stellar surface specify the values of g'(qo) and a through 
Eqs. (|4|l. The value of a depends on the radial and azimuthal 
field components at the surface and must be regarded as fixed. 
Conversely, the value of bo/co is pratically independent of the 
value of g'(qo) for bo/co < - 1 .0. This is illustrated in Fig.|2]and 
is discussed further in Sect. 14.21 In other words, the ratio bo/co 
is only marginally constrained by the boundary conditions and 
can be assumed to vary in order to decrease the total helicity of 
the field. The minimum helicity is obtained for b^/cQ -oo, 
which gives g'iqo) - -8.2133 while for the initial field config- 
uration with /jo/co = -1-1, g'(qo) - -8.2021, coiTesponding to 
a variation of the boundary condition by only 0.13 percent. This 
state of the field has E/Ep = 1.0180 and |//r| = 54.6421 B^R'^, 
i.e., an helicity variation A//r - 0.67B^^ with respect to 
the initial state. Assuming Bq = 10 G and R - 8 x 10^ 
m, the transition from the initial state to this final state re- 
leases an energy AE - 0.0093£'p = 1.6 x 10^'' J. A crucial 
point is the timescale for the release of AE because it deter- 
mines the power available to explain the observed phenomena. 
It depends on the time scale for the dissipation of the helic- 
ity and can be computed from the helicity dissipation rate as 
Thd ~ AH^IidH^ldt) ~ 3.3 x 10^ s, where we made use of 
Eq. (fT2l l. Therefore, the available maximum power, estimated 
from AE/tm, turns out to be ~ 4.9 x 10^° W which is of the 
right order of magnitude to account for the chromospheric hot 
spot and the enhancement of X-ray flux. 

It is important to note that our argument provides us only 
with a lower limit for the energy that can be released because 
we approximate the process as a transition between two lin- 
ear force-free fields. As a matter of fact, when the field is in a 
non-linear force-free state, its energy is greater than the energy 
of the initially unperturbed linear force-free field, allowing the 
system to release more energy than we have estimated above. 
This can explain the modulation of the emitted flux with the 
orbital phase of the planet if only a portion of the stellar corona 
surrounding the flux tube connecting the planet with the star is 
involved in the energy release process at any given time. 

The initial field configuration needs to be restored on a time 
scale comparable with the orbital period of the planet to achieve 
a quasi-stationary situation in agreement with the observations. 
An order of magnitude estimate of the helicity flux coming 
fr om the photospheric motions can be obtained from Eq. (11) 
of iHevvaerts & Priest (Il984) as 



dt 



4jTR^Blv^a'\ 



(13) 



where Vg is the velocity of magnetic flux emergence at the pho- 
tosphere, that we can take as a fraction, say 0.1, of the con- 
vective velocity, yielding Vg - 150 m s"'. The timescale for 
restoring the helicity of the above magnetic configuration is: 



A//r 



Thr = 



3.2 X 10^ s. 



(14) 



dH^/dt 

where dH^ildt comes now from Eq. ( fT3] l. It is comparable to 
the fastest helicity dissipation timescale applied above and is 
significantly shorter (i.e., a; 10 percent) than the orbital period 
of the planet. The magnetic energy flux associated with the he- 
licity build up can be estimated as: ~ InR^B^v^j ^ ~ 4.8 x 10^" 



1 5 

(r/R) 



Fig. 1. The modulus of the magnetic field on the equatorial 
plane of the star normalized to its value at the surface vs. the ra- 
dial distance from the star for bo/co - -1.1 and a = -0.127?"' 
(solid line). For comparison, we plot also the case of a potential 
field (dotted line) and of a field decreasing as r"^ (dashed line). 



W, which can account for the observed energy release in the 
corona. 

It is interesting to note that the greater the magnetic field 
of a hot Jupiter, the more efficient the triggering of helicity 
dissipation in a stellar corona by the planet itself because the 
dissipated power scales as in Eqs. (O and ( fTOl l. However, 
even if the planetary magnetic field is negligible, we still ex- 
pect some dissipation by the currents induced in the planetary 
conductiv e interior by its motion through the stellar coronal 
field (Lai ne et al.ll2008h or by the currents associated with the 
Alfven wave s excited by the motion o f the planet through the 
stellar wind /Preuss e et al.ll200 ?. '2006V The corresponding he- 
licity dissipation rate in Eq. JTOb is expected to be at least two 
orders of magnitude smaller than when the planet has a field Bpi 
of 5 - 10 G, essentially because the cuiTent dissipation is con- 
fined within a much smaller volume. Therefore, the efficiency 
of the proposed mechanism is likely to be significantly reduced 
when the planetary field vanishes. 



4.2. The intermittent nature of SPMI 

The fact that chromospheric hot spots rotating sychronously 
with the planetary orbit have not always been observed in 
HD 179949 and u And (Shkol nik et al.l2008h can be interpreted 
by assuming that their activity cycles are very short, i.e., of 
the order of 1 - 2 years. If the large scale dipole field cycli- 
cally reverses its direction, we have alternate phases of strong 
and weak interac tions with the magnetic field of the planet 
(see Lanzal 2008 ). Such a hypothesis has received some sup- 
port from the observation of two consecutive field reversals in a 
couple of years in t Boo (iDonati et al.,,2008; .Fares et al .„2009h . 
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Fig. 2. The first derivative of the function g at the surface of the 
star (i.e., r = R) vs. bo/co for a = 0.127? 



Moreover, ICranmer & Saail (|2007l) and IShkolnik etaP (l2009h 
build a statistical model for SPMl, based on the observations 
of the solar photospheric field along cycle 22 and its extrapo- 
lation to the corona by means of a potential field model, which 
predicts an intermittent interaction. The relative durations of 
the on and off phases depend on the phase of the solar cycle 
and are in general agreement with the few available observa- 
tions reported in Sect.|2] Nevertheless, another interpretation is 
possible which does not require a global reversal of the mag- 
netic field or a complex geometry of the coronal field, but only 
a change of the photospheric boundary conditions that induces 
a change of the global topology of the coronal field. 

We shall explore the effect of a change of the bo undary 
condit ions on the force-free configuration adopted by iLanza 
(l2008h to model SPMI in the case of HD 179949. It assumes 
a = -0.12i? ' and bo/co - -1.1. The mean of the merid- 
ional field component over the surface of the star B^g ^ in the 
rh.s. of the second of Eqs. (|4]l fixes the value of g'iqo) from 
which we can derive the model parameter bo Icq. A plot of 
g'ilo) vs. bo/cQ is given in Fig.|2]for a - -0.12/?"' and shows 
that there is a large interval of bo/cQ where the variation of 
g'iqo) is very small. In other words, a modest variation of the 
meridional field component produces a remarkable variation of 
bo/ Co- On the other hand, to change significantly the param- 
eter a we need a remarkable change of the axisymmetric az- 
imuthal field component over the surface of the star that can 
be achieved only on timescales comparable with the stellar ac- 
tivity cycle. Therefore, the dependence on the boundary condi- 
tions suggests to explore the modification of the field topology 
for a variable bo/co fielding a fixed. 

The plots of the radial function g(q) and its first deriva- 
tive g'iq) vs. q are shown in Fig.[3]for several values of bo/co 
and a = -0.12R-\ For (bo/co) < -0.9745, the function g is 
monotonously decreasing and all the field fines have both ends 




\ -4 



Fig. 3. Upper panel: The radial function g vs. q for a - 
-0.12/?"' and different values of bo/co as indicated by the 
different linestyles, i.e., solid: bo/co = -0.107; dotted: -0.2; 
dashed: -0.5; dash-dotted: -0.9745 (ciitical value for the for- 
mation of the azimuthal flux rope); dash-dot-dotted: -2.0. 
Lower panel: The first derivative g' vs. q for a - 0.12/?"' and 
different values of bo / co, according to the same linestyle coding 
adopted in the upper panel. 



anchored onto the photosphere. A sketch of t he meridional sec- 
tion of the field lines is given in Fig. 1 of Lanzal (2008) for 
bo/co = -1.1. On the other hand, for (bo/co) > -0.9745, 
the derivative g' vanishes at two points within the interval 
qo < q < qh, and the field develops an azimuthal rope of flux 
centred around the maximum of g{q) and whose inner radius 
coincides with the minimum of g(q). The dependences of the 
magnetic field energy E, outer radius tl, and absolute value 
of the relative helicity |//r| on bo/co for a - -0.12/?"' are 
illustrated in Fig. H] For -2.0 < bo/co < -0.3, the energy, 
outer radius, and absolute value of the relative helicity all in- 
crease very slowly with the parameter, thus we restrict the plot 
to the interval showing the steepest variations. The parameter 
corresponding to the Aly energy is bo/co - -0.107 and it is 
marked by a vertically dashed line. A meridional section of the 
field lines when the field has the Aly energy is plotted in Fig.|5] 
and shows how the azimuthal flux rope has extended to occupy 
most of the available volume, while the domain with field lines 
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Fig. 4. Upper panel: The magnetic field energy as a frac- 
tion of the corresponding potential field energy vs. feo/co for 
a - -0.127? '. The Aly energy limit is marked by the hori- 
zontal dotted line while the vertically dashed line marks the 
corresponding value of i>o/co - -0.107 and has been reported 
into all the other plots. Middle panel: The outer radial limit 
of the field vs. /jq/cq for a - -0.12/?"'. Lower panel: the 
absolute value of the relative magnetic helicity vs. bo Icq for 
a = -Q.nR-\ 

connected to the surface of the star has been squeezed below a 
radial distance ~ 5R. 

We have illustrated the dependences of g' , field energy 
and helicity on the parameter bo Icq for the particular case of 
a - -0.12/?"', but their qualitative behaviours are the same 
also for different values of a. Specifically, we have explored 
numerically the dependence of EIE^ and |//r| in the rectangle 
(0.01/?-' < lal < 0.6/?-') X (-1.5 < boIcQ < -0.001) find- 
ing that they are always monotonously increasing functions of 
bo I Co for any fixed value of a. 

In Sect. 14.1.2! we have assumed that the helicity of the 
stellar corona is determined by a dynamical balance between 
the opposite contributions of the emerging magnetic fields and 
photospheric motions that build it up, and the orbital motion 
of the planet that triggers a continuous dissipation of helic- 
ity and magnetic energy in the corona. Therefore, the photo- 
spheric boundary conditions that fix the value of the helicity in 
our model can be regarded as a result of those processes that 
rule the helicity balance of the stellar corona. 

To explain the transition between states with and with- 
out a chromospheric hot spots, we assume that the helicity of 
the coronal field and the corresponding boundary conditions 
can vary on a timescale shorter than the stellar activity cycle. 
Specifically, the flux of helicity into the corona and the merid- 
ional component of the surface field are expected to vary as 
a result of the fluctuation s characterizing turbulent hydromag - 
netic dynamos (cf., e.g., [Brandenburg & S ubrama nianll2005b . 
As a consequence, boico varies, spanning a certain range of 
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Fig. 5. Meridional section of the magnetic field lines for a - 
-0.12/?-' and boIco - -0.107 which corresponds to the Aly 
energy limit. The distance from the rotation axis is indicated 
by 5 = r sin 6, while the distance from the equatorial plane is 
z = r cos 0, with /? being the radius of the star. 



values. If such a range is large enough, the value of the pa- 
rameter can sometimes cross the thresold boIco = -0.9745 for 
the transition to a flux rope topology, thus halting the SPMI 
mechanism described in Sect. 14.1.21 Conversely, when the pa- 
rameter crosses the threshold in the reverse direction, SPMI 
is resumed. The relative durations of such on and off phases 
depend on the statistical distribution of the fluctuations of the 
meridional component of the surface field which is presently 
unknown. However, an on/off transition can take place on a 
time scale as short as 10^ - lO*" s assuming the helicity fluxes 
and dissipation rates estimated in Sect. 14.1.21 

The phase lag A(p between the planet and the c hromo- 
spher ic hot spot depends on boIco for a fixed a (see iLanza 
20081 for the method to compute A(p). We plot such a de- 
pendence in Fig. |6] When the field has no flux rope, i.e., 
bo I Co < -0.97, A(p depends only slightly on boIco varying only 
by +20° around the mean observed value of ~ 70°, even for 
values of boIco as small as -3.0, i.e., well beyond the range 
considered to model the observations of HD 179949. The vari- 
ation of A(p becomes steep only when the field is very close to 
develop a flux rope. The inner radius of the flux rope reaches 
the distance of the planet (r//? - 7.72) when boIco - -0.908. 
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Fig. 6. The phase lag between the planet and the synchronous 
chromospheric hot spot vs. feo/co for a - 0.12/? 



We assume that the transition from, say, bo/cQ - -1.1 to -0.9 
occurs on a time scale so short that it has no observable conse- 
quences on A(p given the limited duty cycles of present ground- 
based observations. 

By increasing bo/co, the energy of the field grows up to the 
Aly limit where an instability opening the field lines may pos- 
sibly be triggered. The magnetic field intensity vs. the distance 
from the star is plotted in Fig.|7]for a model with a = -0.12/? ' 
and the Aly energy. Following a rapid initial decrease, closely 
similar to that of a potential field, the field stays almost constant 
for 5 < (r/R) < 14, and then decreases slowly far away from 
the star This is due to the slow decrease of the field intensity 
with radial distance inside the flux rope. A first consequence 
is a remarkable increase of the power dissipated at the bound- 
ary of the planetary magnetosphere, as given by Eq. In the 
case of HD 179949, the field intensity at a = 7J2R is ~ 3 times 
greater than in the absence of the rope (cf. Fig.IIJ^ leading to an 
increase of Pa by a factor of ~ 4. Nevertheless, the released en- 
ergy cannot reach the chromosphere and no hot spot is formed 
in the present case. If a balance is eventually reached between 
the helicity flux from the photosphere and the dissipation in the 
corona, the flux rope configuration may become stationary. Its 
end may occur either by an increase of the energy above the 
Aly limit which may open up the field lines, or by a sudden 
decrease of the helicity flux which will change the topology of 
the field into one with all field lines connected to the stellar 
surface, thus resuming a chromospheric hot spot. 

The explanation suggested above for the intermittent be- 
haviour of SPMI can be tested by measuring the direction of 

(s) 

B'g', i.e., the meridional field component at the stellar surface, 
whi ch is possible by means of spectropolar imetric techniques 
(cf. iMoutou et all 120071 lOonati et al1l2008h if the star rotates 
fast enough (y sin / > 12 - 15 km s"^). If the on/ofi" SPMI tran- 
sition is not associated with a reversal of Bg \ the present ex- 




Fig. 7. The modulus of the magnetic field on the equatorial 
plane of the star normalized to its value at the surface vs. 
the radial distance from the star for /jq/cq = -0.107 and 
a - -0.12/? ' (solid line). For comparison, we plot also the 
case of a potential field (dotted line) and of a field decreasing 
as r"^ (dashed line). 



planation gains support. In principle, by measuring the surface 
components of the field and the angle A<p with sufficient ac- 
curacy, it is possible to estimate bo/c^, thus constraining the 
model parameters during the SPMI on phases. Specifically, we 
can use the first and the third of Eqs. Q to find Bq and a, and 
then fix the range of b[)/co that reproduces the observed (cf., 
e.g., Fig.|6l). The second of Eqs. (|4]i can be used as an indepen- 
dent check of the accuracy of the linear force-free assumption 
because the resulting value of g'iqo) should be close to -8 if 
our model is indeed applicable. 

4.3. Azimuthal flux rope and planetary evaporation 

An interesting property of the flux rope topology is the pos- 
sibility of storing matter in the stellar corona. The evapora- 
tion of hot Jupiters under the action of the ionizing radia- 
tion of their host stars originates a flow of cool plasma at a 
temperature of ~ 10"* K that escapes from the planetary at- 



mospheres into the stellar coro nae (cf., e.g., lEhrenreich et al 



20081: iMurray-Clay et alJl2009 ). The flux rope geometry keeps 
the evaporating plasma confined into a torus in the equatorial 
plane of the star, thus r naking its detection easier in the case of 
transiting hot Jupiters ('Vidal-Madi ar et alj|2003h . 

Considering a moderately active star and assuming a life- 
time of ~ 300 days for the flux rope configuration accord- 
ing to the duration of the off" phases of chromospheric inter- 
actio n, the evaporated mass is of the order of 5 x 10'^ kg (cf., 
e.g., iMurray-Clay et al.l2009l) . Such a material may eventually 



condense in the form of prominence-like structures around the 
minimum of the gravitational potential inside the flux rope, that 
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is in the equatorial plane of the star, possibly closer to the star 
than the hot Jupiter (cf. Fig.|5]l. It is unlikely that all the evap- 
orated matter collects into a single condensation. The thermal 
instability of the plasma is expected to lead to many conden- 
sations with typical lengthscales of 10^ - 10^ m, comparable 
in orde r of magnitu de to those observed in, e.g., solar promi- 
nences (lFieldlll965h . Moreover, a high degree of inhomogene- 
ity is expected inside each condensation, in analogy with the 
filamentary structure of solar prominences, because the pres- 
sure scale height at a temperature of ~ 10"* K is only of the 
order of 10^- 10^ m. 

Assuming a mass of IMq for the star and a radius R - IRq, 
the gravitational potential energy of the whole prominence ma- 
terial is ~ 2 X 10^^ J, for a mean radial distance of IR. It is only 
1 percent of the magnetic energy of the coronal field, therefore 
it is not expected to perturb appreciably the field geometry. 

On the other hand, in the case of a very active star, i.e., hav- 
ing an X-ray lumi nosity 2-3 orde rs of magnitude greater than 
the Sun (cf., e.g., IPenz et alJl2008 ). the planet evaporation rate 
may increase by a factor of ~ 100, leading to a potential energy 
of the whole prominence material of the order of 10^^ - 10^^ J. 
This can significantly perturb the field, whose configuration 
can no longer be assumed force-f ree. A non-force- free model 
m ay be a pplicable, such as that of Neukirchl ( 1995 ). discussed 
by iLanza (.2008.) . It gives rise to the same field topologies of 
the present force-free model because it is obtained from a lin- 
ear force-free model by means of a linear transformation of the 
independent variable q of the radial function g(q) (cf. Lanzal 
20081). A large amount of prominence material may help to sta- 
bilize the flux rope configuration, because lifting up the promi- 
nences to open or change the geometry of the magnetic field 
lines requires a comparable additional amount of energy. This 
may imply that the flux rope topology is the most stable in the 
case of very active stars, leading to a lower probability of ob- 
serving a chromospheric hot spot synchronized with the planet. 

Prominence-like structures have indeed been detected 
around some young and highly active dwarf stars, through 
the absorption transients migrating across their Ha, Ca II 
H & K and Mg II h & k emission line profiles (e.g.. 
Collier Cameron & Robinsonl[l989) . Typical masses are in the 
range (2 - 6) x 10^"^ kg for flie young (=a 50 Myr), single KO 
dwarf AB Dor jCollier Cameron et al. fll990'). We expect a re- 
markable increase of the number of such phenomena in highly 
active stars hosting evaporating hot Jupiters. Assuming a life- 
time of ~ 300 days for the flux rope configuration, a total evap- 
orated mass of 2 X 10'^ kg can lead to the formation of ~ 100 
prominence-like structures akin those observed in AB Dor. 

The most favourable geometrical conditions for their de- 
tection are found in transiting systems where the promi- 
nences are expected to transit across the disc of the host star 
The radial distance of a condensation inside a flux rope is 
of several stellar radii, in contrast to a minimum distance 
of a few stellar radii observed in the case of the promi- 
nenc es formed by evaporation from the chromosphere of the 
star (iCoUier Cameron & Robinson 1989). This implies that the 
transit times of the absorption features across the line pro- 
files are shorter by a factor of 2 - 5 in the case of flux rope 
condensations. On the other hand, the expected relative de- 



crease of the intensity along the Ha line profile is between 
5 and 20 percent, with typical equivalent widths of the ab- 
sorption features between 50 and 500 mA, as reported by 
Colher Cameron & RobinsonI ( 11989) for AB Dor, whose coro- 
nal parameters can be considered fairly typical of those of very 
active stars. 

Finally, note that a prominence does not form at the dis- 
tance of the planet when the field topology does not contain 
a flux rope because in that case all the evaporated matter will 
fall onto the star, except when the star rotates so fast that the 
centrifugal force reverses the effective gravity at the top of the 
loops interconnecting the planet with the star. 



4.4. Consequences for X-ray emissions 



In our model, all the magnetic field lines are closed, so the 
corona consists of closed loops. However, if we adopt a more 
realistic model, s uch as the non-linear force-free model of 
Iplyer et al] (l2004l) . it is possible to have open field lines which 
may account for the configuration observed in solar coronal 
holes. Also in those non-linear models we can have a coronal 
field with a large azimuthal flux rope which affects the topol- 
ogy of the coronal field lines close to the star. From a qualitative 
point of view, we can refer to Fig. |5] showing that all the mag- 
netic structures at low and intermediate latitudes must have a 
closed configuration with a top height lower than ~ 4-5 R be- 
cause they must lie below the flux rope. This implies that open 
field configurations, akin solar coronal holes, are not allowed 
at low latitudes when a sizeable flux rope has developed in the 
outer corona. Since closed magnetic configurations are charac- 
terized by X-ray fluxes up to ~ 10-100 times greater than 
coronal holes, an azimuthal flux rope configuration is expected 
to be associated with a greater X-ray luminosity of the star than 
in the case when all magnetic field lines are connected to the 
photosphere. Of course, this is independent of the presence of 
a planet. 



When a star is accompanied by a close-in planet, it may 
increase the dissipation of magnetic h elicity and energy, a s 
conjectured in Sect. 14.1.21 According to lKashyap et al.l (l2008h . 
stars with a distant planet have an average X-ray luminos- 
ity Lx ~ 7 X 10^'^ W, while stars with a hot Jupiter have 
Lx ~ 5.5 X 10^' W. Considering the model introduced in 
Sect. 14.1.21 the dissipated power is proportional to B^, where Bq 
is the photospheric magnetic field (see Sect.O. This holds true 
both for a field topology with all field Unes connected to the 
photosphere and with an azimuthal flux rope. We find that for 
both field topologies an average photospheric field Bq ~ 30 G 
is sufficient to account for the enhancement of X-ray luminos- 
ity in stars with hot Jupiters. Since our model underestimates 
the free energy available in real non-linear force-free fields, 
such a value of Bq should be regarded as an upper limit, thus 
our estimate agrees well with the available observations (e.g., 
iMoutou et al.,.2007i) . 
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4.5. Consequences for radio emission 

The frequency of radio emission from exoplanets depends on 
the planetary magnetic fields (Zarka 2007). According to the 
scaling laws adopted bv lGriefimeier et alj (12004 ). the magnetic 
moments of hot Jupiters are expected to be one order of mag- 
nitude smaller than that of Jupiter owing to tidal synchroniza- 
tion between their rotation and orbital motion. Adopting a sur- 
face field of 10 percent of Jupiter, i.e., 1.4 G, the cyclotron 
maser emission should peak at ~ 4 MHz, which is impossi- 
ble to detect because it falls below the plasma frequency cutoff 
( Jardine & Collier Camerorj2008 ). The situation is much more 
favourable if the magnetic fields of hot Jupiters are indeed as 
strong as predicted by the models of IChristensen et ak I (l2009h 
in which the planetary dynamo is powered by the internal con- 
vective motions and the field intensity is pratically independent 
of the rotation rate of the planet. Assuming a field intensity at 
the surface of 25 G, the cyclotron maser emission peaks at fre- 
quencies around 70 MHz. 

In addition to the emission from the poles of the planet, we 
expect radio emission from the stellar corona. When the field 
lines are connected to the stellar surface, electrons accelerated 
at the reconnection sites travel down to the star producing emis- 
sions up to the GHz range from localized regions above photo- 
spheric spots with fields of 10^ - 10^ G. However, this is not 
possible when the field has a topology with an azimuthal flux 
rope. In this case the coronal field intensity is of the order of 
10"^ G and the emission peaks at very low frequencies, compa- 
rable to or below the plasma frequency, i.e., it is self-absorbed 
before escaping from the emitting region. In this case, we ex- 
pect detectable radio emission only from the poles of the planet, 
i.e., from a highly localized region. The beaming of cyclotron- 
maser emission may additionally decrease the power toward 
the observer thus explaining the lack of detection at frequen- 
cies of 50- 100 MHz. 

In conclusion, we expect that the best chances of detecting 
radio emission from exoplanets with the current instrumenta- 
tion can be achieved when observing systems with a chromo- 
spheric hot spot synchronized with the planet. In this case, a 
fraction of the electrons accelerated at the reconnection sites 
inside the long loop connecting the star with the planet may be 
driven to the polar regions of the planet and to the stellar sur- 
face increasing the irradiated power In this case, the emission 
is expected to be strongly modulated with the orbital period of 
the planet. 



5. Discussion 

We have presented a model for the coronal magnetic field of 
late-type stars that allows us to investigate the magnetic inter- 
action betwe en a sta r and a close-in exoplanet. The same model 
was used bv iLanzal (2008) to explain the observed phase lags 
between chromospheric hot spots attributed to SPMI and the 
planets. Now, we have discussed the energy budget of SPMI 
and its intermittency on the base of that model. We suggest that 
the coronal field evolution in a star hosting a hot Jupiter is ruled 
by a dynamical balance between the helicity coming up into the 
corona from the photosphere and that dissipated by the recon- 



nection events triggered by the orbital motion of the planet. At 
any given time, most of the energy is dissipated in the loop con- 
necting the planet with the stellar surface and in its neighbour 
magnetic structures rather than at the boundary of the planetary 
magnetosphere. 

This scenario is quite different from that characteristic of 
stars with distant planets, like our Sun. In that case the accu- 
mulation of magnetic helicity and energy in the corona leads 
to an instability of closed field structures that erupt as coro- 
nal mass ejections (CMEs) eliminating the excess of helicity 
(IZhang & Lowll2005l: IZhang et al. I l2006h . ■ n the present model, 
in addition to the CME mechanism, the interaction between 
the coronal field and the planetary magnetosphere takes part in 
the helicity dissipation process. In principle, one expects that 
a large flare that extends over most of the stellar corona may 
sometimes occur thanks to the capability of the planet to trigger 
a large-scale helicity dissipation process. If most of the coronal 
helicity is dissipated during such an event, the maximum avail- 
able energy can approach the difference between the Aly limit 
and the energy of the potential field with the same radial com- 
ponent at the photosphere, i.e., AE^^^ = 0.66£'p. For the case of 
HD 179949 with an assumed photospheric field of Bo = 10 G, 
we have AE^ax = 1.1 x 10^^ J. Assuming a photospheric mean 
field Bo = 40 G and R - 0.75 R©, as suggested by the observa- 
tions of the K dwarf HD 189733 bv Moutou et al.. (2007) . we 
get A/imax = 5.2 X 10-^ J. Such large energies may explain the 
superflares observe d in some dwarf stars, g iving support to a 
conjecture by Rubenstein & Schaefer (l2000h . 

Processes like CMEs, i.e., capable of reducing the helic- 
ity of the stellar field, are require d for the operation of a stel- 
lar hydromagnetic dynamo (.Blackman & Brandenburg 2003 



Brandenburg & Subramanian 2005b . Therefore, a hot Jupiter 
may help the star to get rid of the helicity generated by its dy- 
namo in the convection zone, increasing dynamo efficiency and 
the overall level of magnetic activity. This may ultimately ex- 
plain the greater X-ray luminosity of stars with hot Jupiters. 
Moreover, shorter activity cycles may be expe cted (cf., e.g., 
§ 11.2.1 of Brandenburg & SubramanianI 2005 ). which might 



account for the 2-year magnetic c ycle suggested by the latest 
observations of t Boo (IFares et al.ll2009) . 



The modulation of the helicity loss with the orbital period 
of the planet might account for the photospheric cool spots that 
ap pear to rotate synchronously with the planet, as conjectured 
by iLanzal (l2008h . These spots are related to the emergence of 
magnetic flux from the convection zone that may contribute 
to the formation of hot spots in the chromosphere and in the 
corona by reconnecting with pre-existing fields, thus contribut- 
ing to the energy budget of SPMI. 

Our model predicts a field intensity that decreases slower 
than that of a potential field far away from the star. In the case of 
the flux rope topology, the field decreases even slower than r"^ 
in the outer part of the closed corona. This may enhance the dy- 
namical coupling between the star and the planet, as suggested 
by some preliminary computations of the angular momentum 
exchange between CoRoT-4a and its hot Jupiter having an or- 
bital semimajor axis a ~ 17.47?, where R is the stellar radius. 
They indicate that the rotation of the outer convection zone of 
the star may have been synchronized with the orbital motion of 
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the planet if the surface field Bo> 10 G and the age of the sys- 
tem is > 0.5 Gyr. These results may explain the tight synchro- 
nization observed in the CoRoT-4 system which is impossible 
to understand with tidal models given t he large distance of the 



Lanza et alj|2009bh . 



hot Jupiter and its low mass (~ 0.7 Mj, 

A limitation of the present approach is the use of a linear 
force-free model for the stellar coronal field. This mainly af- 
fects our energy estimates, while the main topological features 
of our hnear model are shared by no n-linear models (cf., e.g., 
Fiver et aljr2004l; IZhang et alj|2006l) . The latter can be useful 
to treat the contribution of the outer coronal fields because they 
can provide us with field configurations that extend to the in- 
finity with a finite magnetic energy. However, given the much 
greater mathematical complexity of non-linear force-free mod- 
els, the present treatment is preferable for a first description of 
the most relevant physical effects. 

Non-linear models can be useful also to investigate the 
stability of coronal magnetohydrostatic configurations. In the 
framework of the adopted model, stability is warranted by 
Woltjer theorem because a linear force-free configuration is 
the minimum- energy state f or given total helicity and bound- 
ary conditions (lBergeJl985h . However, a real coronal field that 
is in a non-linear force-free state may become unstable before 
reaching the Aly energy limit by, e.g., kink modes, when it de- 
velops an azimuthal flux rope. 

6. Conclusions 

We ha ve further investigated the model proposed by Lanzal 
(12008') to interpret the observations of star-planet magnetic in- 
teraction. A linear force-free model of the stellar coronal field 
has been applied to address the energy budget of the interac- 
tion and to understand its intermittency. We propose that the 
magnetic helicity budget plays a fundamental role in the inter- 
action. An hot Jupiter contributes to this budget by increasing 
the helicity dissipation that triggers an additional magnetic en- 
ergy release in the stellar corona. 

The transition between phases with and without a chromo- 
spheric hot spot rotating synchronously with the planet is inter- 
preted as a consequence of a topological change of the coronal 
field induced by an accumulation of helicity. Its timescale de- 
pends on the mechanisms ruling the helicity fluctuations in the 
stellar hydromagnetic dynamo and the helicity budget of the 
stellar corona which are presently poorly known. However, a 
timescale as short as 10^ - 10^ s could in principle be possible 
(cf. Sect.gJli. 

The model can be tested in the case of sufficiently rapidly 
rotating stars, such as t Boo or HD 189733, by combining spec- 
tropolarimetric measurements of the stellar photospheric fields 
with Ca II K line observations to determine the phase lag be- 
tween a planet-induced hot spot and the planet. Such simul- 
taneous measurements can in principle constrain the topology 
of the field by allowing us to estimate the parameters of the 
force-free field model as detailed at the end of Sect. 14.21 

The present model also bears interesting consequences for 
the coronal emissions of stars hosting hot Jupiters. It may ex- 
plain why stars with a close-in giant planet have, on the aver- 
age, a higher X-ray luminosity than stars with a distant planet. 



Moreover, it suggests that the best chances to detect radio emis- 
sion from hot Jupiters or their host stars are found in systems 
showing a chromospheric hot spot rotating synchronously with 
the planet. 

We have also investigated the consequences of the different 
field topologies for the confinement and the storage of the mat- 
ter evaporated from a planetary atmosphere under the action of 
the radiation from the host star When the field has an azimuthal 
rope of flux encircling the star, the evaporated matter can nei- 
ther escape nor fall onto the star and is expected to condense in 
the outer corona forming several prominence-like structures. It 
can in principle be detected in the case of rapidly rotating and 
highly active stars through the observations of transient absorp- 
tion features moving across the profile of their chromospheric 
emission lines. 

Simultaneous optical. X-ray and radio observations can 
prove the association between a chromospheric hot spot in- 
duced by a hot Jupiter and the X-ray and radio emission en- 
hancements expected on the basis of our model because the 
energy is mainly released in the corona of the star and then 
conveyed along magnetic field lines to heat the lower chro- 
mospheric layers. Most of the coronal energy should be re- 
leased within one stellar radius, where the magnetic field is 
stronger Therefore, no large phase lags are expected between 
the chromospheric enhancement and the X-ray and radio en- 
hancements varying in phase with the orbital motion of the 
planet. On the other hand, when there is no signature of SMPI 
in the chromosphere, we expect that also the modulation of 
the X-ray flux with the orbital motion of the planet is sig- 
nificantly reduced. No detectable radio emission is expected 
in this case, except when the planet has a polar field of at 
least 20 - 30 G. Finally, we expect to observe the signatures 
of several prominence-like condensations in the coronae of 
rapidly rotating (vsin/ > 40 - 50 km s"'), highly active stars 
(Lx ~ 10^^ - lO^-' W) hosting transiting hot Jupiters, if our 
assumption that coronal flux ropes have a typical lifetime of 
200 - 300 days is true. 
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